Causal structure of the entanglement renormalization ansatz 
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We show that the multi-scale entanglement renormalization ansatz (MERA) can be reformulated 
in terms of a causality constraint on discrete quantum dynamics. This causal structure is that of 
de Sitter space with a flat spacelike boundary, where the volume of a spacetime region corresponds 
to the number of variational parameters it contains. This result clarifies the nature of the ansatz, 
and suggests a generalization to quantum field theory. It also constitutes an independent justifica- 
tion of the connection between MERA and hyperbolic geometry which was proposed as a concrete 
implementation of the AdS-CFT correspondence. 



The multi-scale entanglement renormalization ansatz 
(MERA), introduced in has been successfully used to 
model the physics of many low-dimensional strongly cor- 
related quantum many-body systems [H-Q, thanks to its 
ability to simultaneously represent correlations at widely 
different length scales. It is defined as the set of states 
which can be output by a quantum circuit with fixed in- 
put and a given gate structure. The components of the 
gates arc the variational parameters. 

In this paper, we show that only one general property 
of the circuit is necessary to define the ansatz: its causal 
structure, which is defined as a constraint on the flow of 
information. By clarifying how and when the specific de- 
tails of the circuit are unimportant, this approach unifies 
and simplifies the formulations of a MERA, and opens 
the way to a deeper theoretical study of the ansatz, such 
as its relation to other theories of renormalization [lol |. 
We show that the causal structure of MERA is essentially 
identical to that of dc Sitter spacetime, i.e. an expanding 
universe undergoing an exponential multiplication of its 
local degrees of freedom, such as the infiationary period 
in cosmology. 

This naturally points to quantum field theory on de 
Sitter space as a continuous generalization of MERA able 
to represent the state of critical quantum field theories. 
We show that this proposal is compatible with that of 
Ref. |ll| . but better constrained, and amenable to the 



tools developed in the context of quantum cosmology. 
Additionally, this result constitutes a connection between 
MERA and hyperbolic geometry (the Euclican form of 
both anti dc Sitter and dc Sitter spacetimes), which is 
independent from the arguments previously proposed in 
connection with the AdS-CFT correspondence [l2l. [l3l|. 



MERA FROM CAUSAL ORDER 

One step of a circuit defining a MERA takes a quan- 
tum state defined on a coarse-grained lattice and isomet- 
rically maps it into the larger Hilbert space of a finer 
lattice. These isometric steps are also required to act 
locally. This implies a finite speed of information prop- 
agation in the circuit. Combined with the exponential 



nature of the successive coarse-graining operations, this 
means that the expectation value of a local observable 
can be evaluated in a time logarithmic in the lattice size. 
In order to see this, note that an expectation value can 
be evaluated by evolving the observable "back in time" 
in the Heisenberg picture and then computing the ex- 
pectation value between the initial fiducial state and the 
resulting observable. When talking about locality it is en- 
lightening to adopt the Heisenberg picture because there 
exists an unambiguous concept of a local observable: one 
which acts nontrivially only on certain lattice sites. The 
locality of the isometrics then implies that an observable 
local to a region S of the lattice is pulled back to an ob- 
servable which is itself local with respect to a region E' 
that is in the causal past of S Q . This map is the dual of 
a local quantum channel mapping states defined on S' to 
states on S. In performing this operation, the rest of the 
isometry can be completely ignored. Furthermore, the 
coarse-graining is such that the Hilbert space dimension 
associated with the causal past S' is no larger than that 
of E, and hence the computational load can only decrease 
at each step, and is independent of the lattice size. 

This suggests that, in defining the ansatz, we could 
replace the ad-hoc requirement that each isometric step 
have a particular gate structure, and just require it to 
pull back local obscrvablcs (on E) to local observables 
(on E'). This property is precisely one of causality as it 
is equivalent to stating that the degrees of freedom out- 
side E' cannot infiuence those inside E through one step 
of the dynamics [l4| . But is such a property sufficient to 
obtain an efficient local parameterization of the isome- 
trics? It was shown by Arrighi et al. [lij in the context of 
a unitary dynamical step U that such causal constraints 
arc sufficient and necessary for U to be implemcntable 
as a circuit of local operations, with some commutativity 
constraints between them. For a MERA, however, we 
need each step to be isometric, for which the causality 
constraints are not sufficient to render the parameteriza- 
tion efficient. Indeed, the isometry can always acausally 
produce a state correlated over arbitrary distances. Here 
we introduce a stronger causality requirement which ap- 
plies to isometric maps as well as to generic quantum 
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channels. We demand that the map can be implemented 
by creating a separable state on extra systems at each 
lattice sites (with no constraint on their dimensions), fol- 
lowed by a causal unitary map. We will say that such 
a map is purely causal rather than just causal. The re- 
sult of Ref. [15[ then trivially extends to this stronger 
requirement, showing that a purely causal channel can 
be implemented with local operations. 

We will represent a given causal relation between two 
lattices as a bipartite graph, where the presence of an 
edge between a site of the input lattice and a site of 
the output lattice indicates a possible causal influence. 
The causality constraints therefore are embodied in those 
pairs of vertices which are not connected by an edge. We 
also interpret a particular causality relation as represent- 
ing the set of isometrics which arc purely causal with 
respect to it. We will show that 
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where time flows upward, and the right-hand side rep- 
resent the set of isometries which can be implemented 
by the circuit (or tensor network) obtained by replacing 
each box by an isometry. This is precisely one step of the 
binary MERA [i|,[3], except for the fact that in our case 
there is no constraint on the dimension of the Hilbcrt 
spaces associated with the intermediate wires. However, 
the fact that each box must be an isometry limits their 
input of dimension to that of their output. 

Below we also show that the ternary MERA Q is 
equivalent to such a natural causality constraint. More 
generally, our prescription together with the construc- 
tive localizability result introduced in Ref. [3l allows for 
the construction of circuits with equivalent properties on 
arbitrary lattices, including lattices embedded in higher 
dimensional spaces. We defer the explanation of how a 
general circuit equivalent to a given causality structure 
can be designed to Ref. fisf . 

Our approach works just as well if we allow each step to 
be implemented by a quantum channel rather than just 
an isometry, hence allowing in principle for the character- 
ization of mixed states with long range correlations, such 
as critical thermal states. If the state to be described is 
classical one may furthermore constrain the local quan- 
tum channels to be stochastic maps (i.e. to map diagonal 
matrices to diagonal matrices). 



CONNECTION WITH DE SITTER SPACE 

Let us call an "event" a lattice site at a given coarse- 
jraining step. Each event is associated with the Hilbcrt 




FIG. 1: Partial ordered set corresponding to the binary 
MERA in one dimension, embedded in such that the speed 
of light is equal to 1 everywhere. The output lattice is the 
top row of dots and time flows upward. The black circles are 
events and the lines segments are causal links. 



space dimension of the corresponding lattice site. The 
causal relations between successive coarse-graining gen- 
erates a partial order between any two events, i.e., A < B 
if A is in the causal past of B. Furthermore the original 
causal relation can be recovered uniquely from the par- 
tial order by noting that it is defined by the causal links: 
pairs of related events with no events "in between", i.e., 
{A, B) forms a link \i A < B and there is no C such that 
A<C <B. 

Therefore, one can recover the MERA simply from the 
causal order between the set of events, together with the 
dimensions of their assigned Hilbert spaces. For instance, 
the usual binary MERA Q is implied by the partial order 
shown in Figure [H Such "causal sets" have been studied 
before as discrete models of spacetime 3, 17 1. The idea 
follows from the fact that the geometry of a manifold with 
Lorentzian signature can be recovered exactly from the 
partial order between events induced by the metric, to- 
gether with the volume form. The causal order directly 
makes sense for a discrete spacetime. For the volume 
form, a natural postulate is that it corresponds to the 
counting of events. In our case, assuming for simplic- 
ity that all events are associated with the same Hilbert 
space dimension, the number of events in a given space- 
time region is proportional to the number of variational 
parameters, thank to the local representability result. 

In order to see what metric a MERA on a d- 
dimensional lattice may correspond to, the easiest is to 
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first parameterize its events by coordinates in which the 
speed of light is constant (and equal to 1), i.e. in a space- 
time with metric ds^ = f{t, xi, . . . , Xd){~dt^ + J2i dxf). 
We suppose that each coarse-graining increases the lat- 
tice spacing by a factor a, and that sites at the (fc -I- l)th 
coarse-graining step have a causal influence on the sites 
of the fcth step within a radius ra^ . Then a constant 
speed of light (equal to 1) is achieved by embedding the 
fcth coarse-graining at time t = —ra^j (a — 1). We choose 
it negative so that it increases chronologically with the 
quantum computation, outputting the final state at time 
to = —f/{a — 1). In order to determine the conformal 
factor /(t, xi, . . . , Xd), wc postulate that in coordinates 
where our lattices are equally spaced in time, and renor- 
malized, the volume form should be constant. This makes 
precise the idea that the number of events in a given re- 
gion of spacetime should be proportional to the volume 
of that region. Such coordinates must be of the form 
T = ^a\og[t/to] and Ci = —fixi/t. The constraint is 
then satisfied by picking f{t) = [a/t)"^. Also, choosing 
to = —r/{a — 1) puts the output boundary fc = at 
r = 0, and /3 = a normalizes the volume element. In the 
coordinates (r, (^i, . . . , Cd) the metric is then 



dT^ ^l^dpdrAr^ dQ. 
a. ^—^ 



where = Q , and the volume form has component 
■\/| det 5I = 1. In the conformally flat coordinates this is 
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FIG. 2: The shaded area is the causal past of two disconnected 
regions of the r = spacelike surface in the static coordinates 
(r, ("). The dashed lines indicate the horizon at |("| = a = 1. 



This analysis yields a possible approach to building 
a continuous MERA: by parameterizing quantum field 
theories on de Sitter space. In order to compare this 
to the "cMERA" proposed by Haegeman et al. [lH, we 
consider a quantum scalar field, with Lagrangian density 

C = i^d^[-g'^>,,,</),. - 

Given our prescriptions, the coordinate system using in 
Ref. [m must correspond to our static coordinate sys- 
tem (r, Ci , • ■ • ) given that the same cutoff is used at each 
time step. This Lagrangian can be quantized by stan- 
dard methods. If 7r(Ci, . . . , C^) are the canonical conju- 
gates to the field operators 0(Ci, ■ • ■ , Cd), then we obtain 
the Hamiltonian H = K + L, where 



Tliis metric is that of de Sitter space. Another common 
coordinate system is given by the time coordinate t to- 
gether with = —axi/to, so that 



ds^ 



-dr 
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Basic properties of the MERA can be deduced 
from considering past lightcones in the coordinates 
(r, Ci, . . . , (^d) with constant volume form. The lightlike 
worldlines can be deduced by applying the coordinate 
change to the Minkowski ones. They are all of the form 

C;(t) = -aM,(l-e^/")-^C.(0)e^/" 

where Ui a unit vector. We see that the causal past of any 
bounded region converges in the infinite past r — ^ — 00 
to the ball of radius a (the cosmological horizon) , which 
contains a fixed number of lattice sites at any given time. 
Figure [2] illustrates this phenomenon. It shows the causal 
past associated with the computation of the correlation 
function between two local observables for d = 1. The 
log-dimension of the tensors that one needs to contract 
at each time step is proportional to the spacclike volume 
of the shaded area at that time. 



K = - j d''c[^' + J2{l,f + v{4>)_ 

i 

is the Hamiltonian for the field on flat spacetime and 

i 

generates the expansion of space. One can easily check 
that, with canonical commutation relations between tt 
and 0, this yields the right Heisenberg equations of mo- 
tion. This structure for H is compatible with that pro- 
posed by Haegeman et al., although our approach yields 
stronger constraints given that it requires exact causality. 

We note also that the Euclidean form of de Sitter space 
is identical to that of anti de Sitter space, namely hy- 
perbolic space. There has been great interest recently 
in MERA as providing a possibly concrete form of the 
AdS-CFT correspondence for its ability to represent crit- 
ical scale-invariant states and the similarity in which en- 
tanglement entro py is calculated using a minimal sur- 
face in the bulk [12 
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This suggests an interpreta- 
tion where the MERA circuit represents a field theory 
on AdS spacetime, whereas the state it describes is a 
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CFT on its boundary. The Euclidean form of our field 
theory lives precisely on the same spacetime, and with 
the same boundary, as the Euclidean form of AdS space- 
time. Therefore our approach provides an independent 
way of relating MERA to hyperbolic geometry. 



is an isometry as can be checked by tracing out A' and 
B' on both sides of Equation ([3]). Then the Stinespring 
dilation theorem tells us that there is an isometry (here 
just a ket) \ip) embedding C into the Hilbert space of the 
system A' such that 



CAUSALITY AND LOCALITY 

We will now sketch the proof of the statement rep- 
resented by Equ. [TJ Namely, that the set of isometrics 
which are purely causal with respect to the causality re- 
lation represented on the left-hand side is equal to the 
set of isometrics which can be formed by replacing each 
box of the right hand side by an isometry (without con- 
straint on the dimension of the Hilbert spaces associated 
with the middle wires). First we consider a unitary map 
U whose inputs are grouped into systems A and and 
outputs are grouped into systems A' and _B',with the con- 
straint that B cannot influence A'. This means that in 
the Heisenberg picture, any operator X acting on system 
A' is mapped to an operator Y acting only on system A, 
i.e. U^XA' 'i-B')U = Ya (8) 1b, which can be rewritten 
as 

{XA'®lB')U = UiYA®lB)- (2) 

This implies that for any pure state of B, 

y = (1 ® {x\) iX(g,i)u{i(g> \x)). 

The trick, inspired by Ref. [l^, is to replace the local 
operator AT by a swap between A' and a new system 
C in Equation If we initialize the system C to an 
arbitrary state \y) and trace it out after the action of U 
and the swap, the left hand side becomes simply U. This 
yields the expansion 
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where the states |a;) and \y) can be chosen arbitrarily. 
This is the only algebraic property that we will need. 
The vertical bar ending the fourth wire means that this 
system is traced out: hence both sides of this equation 
represent channels rather than just operators. The chan- 
nel on the left-hand side is just p ^ U pU'^: being unitary, 
it is a minimal Stinespring dilation of the channel on the 
right hand side. From the unicity of the Stinespring di- 
lation of a channel, the right-hand side has also only one 
Kraus operator. To find its precise form, first note that 
the operator 
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It follows that we can replace the channel J\f{p) := 
tTA'UpW in Equ. © by p X pX'' , with X := 
{1b' <8) {iPIa'W. Furthermore, since the whole expression 
must be unitary, and hence trace-preserving, the opera- 
tor X is isometric when restricted to its possible inputs in 
the circuit, and can therefore be replaced by an isometry. 

This can be used to parameterize the classes of uni- 
tary maps causal with respect to a relation like that of 
Equ. ([1]) as follows: we start by grouping all the inputs 
(resp. outputs) which have the set of children (resp. par- 
ents) to obtain a new causal relation on the grouped sys- 
tems. If the resulting graph is such that removing one 
particular input A breaks it into two independent parts, 
then the remaining inputs and outputs can be grouped 
so as to satisfy the causality relation 

M 

A 

This represents two causality constraints (i.e. missing 
links). By applying the instance of Equ. [3] allowed by one 
of the constraint, and then again on the first instance of 
U in the circuit for the other constraint, we obtain that 
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for some isometry V. This scheme can be applied re- 
cursively on the remaining copies of U, until the circuit 
respects all the causality constraints. 

If we lift the restriction that our computational step 
be unitary, and assume instead that it is an isometry (as 
is required for a MERA) or more generally a quantum 
channel, then we demand that it can be represented by a 
unitary interaction with a local environment, such that 
the unitary map respects the same causal relation. We 
also require that the environment's initial state is sepa- 
rable. We can then apply our procedure to this unitary 
map to show that it has a local representation. In this 
way, one obtains the result express in Equ. dH). This 
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method also works for the ternary MERA, showing that 
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with the same disclaimer about the fact that the dimen- 
sionality of intermediate wires are not constrained, but 
limited by the fact that the boxes must represent isome- 
tries. 

As mentioned in the introduction, for more general 
causality relations, in particular as applied to higher- 
dimensional lattices, one must use the general prescrip- 
tions introduced in Ref. For completeness, we de- 
scribe their construction here. A unitary map respecting 
a general set of causality constraints can be implemented 
locally as follows: one first produce ancillas matching the 
number and dimensions of all the output systems, initial- 
ized with an arbitrary product state. We may then apply, 
for each possible output system, a unitary interaction be- 
tween the ancilla corresponding to that output and its 
parent inputs, i.e. the input from which it is allowed to 
receive information from according to the causality con- 
straints. The only condition on these local unitary inter- 
actions is that they all commute with each other. One 
must then trace out all the input systems, and be left 
with the ancilla, which we identify as the output systems. 
The global unitary respects the causality constraints if 
and only if it can be implemented in this manner. 



ACKNOWLEDGMENT 



[1] 
[2] 

[3; 
[4: 

[5 

[r 

[8 

[9: 

[lo: 
[11 

[12 
[13 

[14 

[15 

[16 

[17 



G. Vidal, Physical Review Letters 99, 220405 (2007). 
G. Evenbly and G. Vidal, Phys. Rev. B 79, 144108 
(2009). 

V. Giovannetti, S. Montangero, M. Rizzi, and R. Fazio, 
Phys. Rev. A 79, 052314 (2009), 0901.4424. 
G. Vidal, ArXiv e-prints (2009), 0912.1651. 
L. Cincio, J. Dziarmaga, and M. M. Rams, Physical Re- 
view Letters 100, 240603 (2008), 0710.3829. 
G. Evenbly and G. Vidal, Physical Review Letters 102, 
180406 (2009), 0811.0879. 

G. Evenbly and G. Vidal, Physical Review Letters 104, 
187203 (2010), 0904.3383. 

P. Corboz, G. Evenbly, F. Verstraete, and G. Vidal, Phys. 
Rev. A 81, 010303 (2010), 0904.4151. 
C. Pineda, T. Barthel, and J. Eisert, Phys. Rev. A 81, 
050303 (2010), 0905.0669. 

C. Beny and T. J. Osborne, in preparation. 

J. Haegeman, T. Osborne, H. Verschelde, and F. Ver- 
straete, Arxiv preprint arXiv:1102.5524 (2011). 
B. Swingle, ArXiv e-prints (2009), 0905.1317. 
G. Evenbly and G. Vidal, Journal of Statistical Physics 
145, 891 (2011), 1106.1082. 

D. Beckman, D. Gottesman, M. A. Nielsen, and 
J. Preskill, Physical Review A 64, 052309 (2001). 

P. Arrighi, V. Nesme, and R. Werner, Journal of Com- 
puter and System Sciences (2010). 

L. Bombelli, J. Lee, D. Meyer, and R. D. Sorkin, Physical 
review letters 59, 521 (1987). 

E. Hawkins, F. Markopoulou, and H. Sahlmann, Classical 
and Quantum Gravity 20, 3839 (2003). 



The author is grateful to Tobias Osborne and Guifre 
Vidal for discussions about this work. This work was 



